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'•‘We  consider  a head-on  collision  between  two  solitary  waves  on 
the  surface  of  an  inviscid,  homogeneous  fluid.  A perturbation 
method  wnich  in  principle  can  generate  an  asymptotic  series  of 
all  orders,  is  used  to  calculate  the  effects  of  the  collision. 

Vie  find  that  the  waves  emerging  from  (i.e.  long  after)  the 
collision  preserve  their  original  identities  to  the  third  order 
of  accuracy  we  have  calculated.  However  a collision  does  ldave 
imprints  on  the  colliding  waves  with  phase  shifts  and  shedding  — 
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of  secondary  waves.  The  latter  propagate  with  diminishing 
amplitudes  in  the  direction  opposite  to  that  of  the  main 
waves.  We  have  also  calculated  the  maximum  run-up  amplitude 
of  two  colliding  waves.  The  result  checks  with  exisiting 
experiments. ^ 
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Abstract 

We  consider  a head-on  collision  between  two  solitary  waves 
on  the  surface  of  an  inviscid,  homogeneous  fluid.  A perturbation 
method  which  in  principle  can  generate  an  asymptotic  series  of 
all  orders,  is  used  to  calculate  the  effects  of  the  collision. 

We  find  that  the  waves  emerging  from  (i.e.  long  after)  the 
collision  preserve  their  original  identities  to  the  third  order 
of  accuracy  we  have  calculated.  However  a collision  does  leave 
imprints  on  the  colliding  waves  with  phase  shifts  and  shedding  of 
secondary  waves.  The  latter  propagate  with  diminishing  amplitudes  in 
the  direction  opposite  to  that  of  the  main  waves.  We  have  also  calcu- 
lated the  maximum  run-up  amplitude  of  two  colliding  waves.  The 
result  checks  with  exisiting  experiments. 
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I.  Introduction 


We  recast  in  Sec.  II,  the  equations  of  motion  of  an  inviscid, 
constant  density  fluid  with  a free  surface  into  a pair  of  equations 
in  terms  of  the  free  surface  elevation  h(t,x,y)  and  the  velocity 
along  a horizontal  stream  bed  W(t,x,y).  These  equations  are  con- 
venient for  a perturbation  scheme  to  be  introduced  in  Sec.  Ill  to 
study  a head-on  collision  between  two  solitary  waves  which  are  small 

in  amplitude  (^—  <<  1),  and  long  in  wave-length  (^—  >>  1).  Where 
o o 

a is  the  measure  of  amplitude,  X the  wave-length,  and  hQ  the  un- 
disturbed depth  of  the  fluid.  The  amplitude  and  the  wave-length 

parameters  are  related  by  Ursell's  ordering  for  theory  of  shallow 
2 3 

water  i.e.  aX  ~ h . We  have  carried  out  the  calculations  to  the 

o 

third  order  of  approximation. 

In  the  first  order  approximation,  we  have  two  independently 

2 v 

moving  solitary  waves  aS(£)  and  bS(n).  S (x)  = sech  -j,  a pro- 
gressive wave  of  permanent  type  which  satisfies  the  Korteweg- 
deVries  equation.  The  variables  £ and  n denote  the  right-  and 
left-going  wave-framed  coordinates  respectively.  The  constants 
a and  b specify  the  heights  of  the  waves.  In  the  second  order 
approximation,  we  find  that  the  wave  field  is  modified  by: 

(1)  quadratic  terms  in  S,  (2)  change  of  wave  speeds  and  (3)  addi- 
tion of  phase  functions  8^(n)  and  <J>^(5)  to  £ and  n respectively. 
For  two  waves  at  large  separation,  i.e.  before  or  after  collision, 
the  first  two  corrections  above  reduce  to  Laitone's  (1960)  second 
order  calculation  of  a single  solitary  wave.  On  the  other  hand, 
we  can  also  use  the  result  in  (1)  to  calculate  the  maximum  run-up 
amplitude  during  the  collision.  Our  result  checks  with  an  earlier 
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calculation  of  Byatt-Smith  (1971).  Since,  up  to  this  order  of 
accuracy,  the  phase  function  0 (or  <j>)  is  a function  of  ri  (or  £) 
alone,  and  differs  by  a constant  value  through  a collision,  this 
value  represents  a phase  shift  for  a right-  (or  left-)  going 
solitary  wave.  This  result  agrees  with  that  of  Oikawa  and  Yajima 
(1973) . In  the  third  order  approximation,  we  obtain  as  before, 
those  three  corrections  to  the  wave  field.  The  first  correction 
is  a polynomial  of  third  degree  in  S.  This  polynomial  in  S and 
the  wave  speed  correction  reduce  as  the  separation  between  waves 
increases,  to  the  known  result  of  Grimshaw  (1971),  Fenton  (1972)  per- 
taining to  a single  solitary  wave  with  third  order  accuracy.  How- 
ever, the  phase  functions  obtained  become  functions  both  of  £ and  n • 
This  gives  different  values  of  phase  shifts  at  different  points  in 
a wave.  Unlike  the  uniform  phase  shift  in  the  second  order  approxi- 
mation which  preserves  the  wave  form,  here  these  variable  phase 
shifts  cause  a distortion  of  the  wave  form  in  addition  to  a simple 
translation.  Since  this  distorted  wave  field  does  not  satisfy  the 
equation  for  waves  which  are  propagated  without  a change  in  shape 
and  speed,  we  are  forced  to  investigate  in  Section  V,  the  slow  time 
evolution  of  this  distorted  wave  forms.  We  find  that  the  waves 
eventually  transform  back  to  their  original  forms  with  shedding  of 
secondarv  wavelets.  These  secondary  wave*  propagate  with  diminish- 
ing amplitude  in  the  direction  opposite  to  that  of  the  main 
waves.  Thus  we  can  still  speak  of  preservation  of  wave  form  up  to 
the  third  order  of  accuracy,  with  the  understanding  that  there  exists 
uniform  phase  shifts  as  well  as  shedding  of  secondary  waves.  We  stop 
our  calculation  at  the  third  order.  However  we  feel  after  working 
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long  with  the  perturbation  scheme  used  here,  that  the  higher  order 
results  will  be  essentially  of  similar  nature.  There  are  two 
imprints  of  a collision  between  two  solitary  waves:  uniform  phase 
shifts  and  shedding  of  secondary  waves.  We  have  also  calculated 
the  maximum  run-up  amplitude  up  to  the  third  order  of  accuracy 
for  two  colliding  solitary  waves.  The  result  agrees  very  closely 
with  the  numerical  and  experimental  results  of  Chan  and  Street 
(1970).  It  also  checks  with  Maxworthy's  experiment  (1976)  for 
a solitary  wave  reflected  from  a wall  (see  Fig.  2). 
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II.  Basic  Equations 

We  take  the  velocity  field  of  our  flow  problem  to  be  described 
by  a potential  4>  (t,x,y , z)  , satisfying  Laplace  equation.  The  free 
surface,  where  the  pressure  vanishes,  is  specified  by  an  unknown 
function  z = h(t,x,y).  The  fluid  is  supported  by  a horizontal 
plane  at  z = 0,  where  we  set  the  normal  velocity  =0.  It  is 

easy  to  show  in  this  case  that  the  potential  <f>  can  be  expressed  as  a 

2 D *!• 

Taylor  series  at  z = 0.  Using  V 4>  = 0 for  z > 0 and  ^ = 0 at  z = 0, 

d Z 

we  obtain  00 


00  2n  „ 

<J>(t,x,y,z)  = l (“Dn  -r^rrr  V^n4> 


n=0 


(2n)  ! 


(1) 


where 


4>(t,x,y)  = <J>  (t  ,x,y,  z=0)  and  V = i 


. 3 


3 x 


+ j , 

3y 


representing  a gradient  operator  in  x,y-plane.  We  can  now  express 
the  kinematic  condition  at  the  free  surface  as 


3 h _ , 3 <j> 

at  + v*'vh  = at  ' 


(2) 


and  the  Bernoulli  theorem  also  applied  at  the  free  surface  as 


+ gh  + j(4>2  + <t>y  + 4>2)  = 


See  Whitham,  Chapter  (13) . 
In  terms  of  4>  these  become 


(3) 


lx  + V* [ J (-l)n  h 
n=0 


2n+l 


(2n+l) ! 


V2n(V<f)]  = 0 , 


(4) 


. gh  ♦ i<w2  + I (-1)"  t 

n=l 

+ \ T V2n"m_14)]  =C  (t) 

m=  0 


(5) 
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where  the  symbol  stands  for  the  vector  inner  product  for  even 

m and  the  usual  arithmetic  multiplication  for  odd  m,  and  (2n)  = 


(2n)  i 


the  binomial  coefficients.  We  shall  consider  waves  in 


m! (2n-m) ! 

a channel.  We  therefore  drop  the  y-dependence  in  h and  . By 
taking  the  x-derivative  of  (5) , we  obtain  in  place  of  (4)  and  (5) 
the  following  two  equations 


3h  , 3 ru  . r , . n h " 3 "w  i n 

at  3 x {hw  £ (-1)  (2n+l)  ! a 2n  ^ 0 

n=l  3x 


, 2n+l  -2n 

h a 


(6) 


- » 2 °°  ,2n  ~ 2n 

3 W , o , W r / i h ro  W 

3t  2 ii'-11  (2n)![3Ux2n-l 


2m 


| T (-Dm(f ) 


m=0 


m , 2n. amw  3 2n  mw 
m axm  ax2n-m 


]}  = o 


(7) 


3 $ 

where  w = ^ represents  the  velocity  at  the  bottom  of  the  channel. 

Both  of  the  above  equations  are  in  the  conservative  form.  It 
is  interesting  to  point  out  that  two  more  conservation  laws  can  be 
derived  from  (6)  and  (7) . They  represent  conservation  of  horizontal 
momentum  and  total  energy,  i.e.. 


,T-r(f  u dx)  + [/  (u2  + ^)dz]  = 0, 

"t  0 x 0 p 


(8) 


gh2+  f (u2+v2)dz]  + / dz  u[-|(u2+v2)  + p + gz]«  (9) 


and 


n ^ ^ ^ 

p = 9(h-z)  + f (g£-  + u + v ^yJv  dz, 
z 


where  u = <f>x  and  v = <Pz  both  are  expressed  in  w by  (1)  . We 


shall  however  in  this  paper  use  Eqs.  (6)  and  (7)  only. 
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III-A.  Perturbation  Solution  for  Two  Head-On  Colliding  Solitary 
Waves . 


Consider  two  solitary  waves,  far  apart,  of  small  but  finite 
amplitude  and  heading  towards  each  other.  The  time  evolution  of 
their  interaction,  and  the  final  state  after  their  collision 
will  be  our  main  concern. 

We  introduce  the  following  coordinates  transformations  (wave- 
frames) 


k(x-CRt) 
n = /z  i (x+C  t) 

O Lj 


(10) 


where  0 < e <<  1,  e is  a dimensionless  parameter  representing  the 
order  of  magnitude  of  the  wave  amplitude.  The  scaling  of  the  hori- 
zontal wave  length  in  accordance  with  Ursell's  relationship  is 
taken  as  /F,  leaving  k and  £ as  the  wave  numbers  of  order  unity 
for  the  right-  and  left-going  waves  respectively.  The  right-  and 
left-going  wave  speed  C , C are  to  be  related  to  the  amplitudes 
of  the  waves.  In  the  limit  of  infinitesimal  amplitude,  they  take 
the  value  of  linear  wave  speed  (ghQ)  ' . Anticipating  a difficulty 

to  show  up  in  our  perturnation  method  we  introduce  the  following 
transformations  of  wave-framed  coordinates  with  phase  functions.* 

* This  is  essentially  the  method  of  strained  cooridnates  first 
introduced  by  Poincare  (1892)  for  ordinary  differential  equations  and 
later  generalized  to  hyperbolic  partial  differential  equations  by 
Lighthi 11  (1949)  and  Lin  (1954).  For  a detailed  discussion  we  refer 
the  reader  to  Van  Dyke's  book  on  Perturbation  methods  in  Fluid 

Mechanics,  chapter  VI.  Academic  Press,  New  York  (1964). 
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SQ  = S - ek€(S,n) 


nQ  = n - e £4>  (S  ,n) 


when  0(S,n)  and  4>  ( C ^ n ) are  to  be  determined  in  the  process  of  our 
perturbational  solution  of  (6)  and  (7) . These  functions,  intro- 
duced for  purpose  or  making  asymptotic  approximation,  allows  us 
to  calculate  phase  changes  due  to  collision. 

Using  (10)  and  (11)  we  obtain  the  transformation  between  de- 
rivatives as 


It  + CR  = K£(CD+CT ) u !_  + ek*  (!&  If  - HItt)  ] 


'R  L' 


an  as  as  an' 


§ _ q ^ = _ ZL(c  +c  ) [k  - + ek£  (—  - - — - — ) 1 

at  k a x d clmk  as  K (as  an  an  as  J 

where  D = (1  - ek£|^)  (1  - e£  f£)  - e2k£  |1  ||-  . 

Setting  h = hQ(l+^),  we  rewrite  (6)  and  (7)  into 

V 

'ft  ±CR,L  + h Ft  ‘ °' 

where  C = l/ghG  is  the  linear  wave  speed,  and 

2 

F+  = +(C-CR  L)  (wfCC)  + j - ± CC w 


« h2n(l+C)2n  82n 

+ l <-Dn  -2 1 1 


. c(i+o  a w 

+ — ~ + 


(2n)  i atax 


2n-l  - 2n+ 


v"1,,  ~ 2n-m. 


1 f , , v m , 2nx  a w 3“ w 

2 ni0  m axm  3x2n_m 


It  is  convenient  to  make  the  following  change  of  dependent  variables 


I 


w + Cc  = 2eCci 
w - C£  = -2eCB , 
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$ 


which  give  us  w = ec(a~3)  and  £ = e(a+3).  Together  with  (12)  and 
(13)  we  can  now  represent  (14)  as 


2Ce(CR+CLH<l  |2  + CMI9,  If11  + (k  lc  + 11  ^ + 


+ - |T(e-*)|^))F+  - o. 


(17) 


A similar  equation  for  3 is  obtained  by  replacing  a by  3 , 

5 by  n , k by  i , and  0 by  <J> . We  need  hereafter  consider  (17)  only, 

We  now  express  the  new  variables  in  the  following  power  series 

2 

a(5,n)  = aQ  + ea^  + e a2  + ... 

3(5/ n)  = 3Q  + e 3 + e2B2  + ••• 


e(5»n)  = eo(n)  + ee1(5,n)  + 


(5»n)  = <(>0  (C)  + e<j>1(5,n)  + ... 

2 2 3 3 

= C(1  + eaR^  +eaR2+eaR3+ 
= C(1  + ebl^  + e2b2L2  + e3b3L3'+ 


(18) 


Substituting  (18)  into  (17)  we  obtain  a lengthy  expression 

in  power  series  of  e,  (See  Appendix).  The  coefficients  of 
2 3 4 

e,  e / e and  part  of  e will,  however,  be  presented  in 
sequel  as  follows: 

i)  Coefficients  of  e. 

The  equations  are  simply 


9a 
( 

9n 


= o 


and 


95 


= 0 


I 


Their  solutions  are 

aQ  = af  ( C)  , BQ  = bg  (n)  (19) 

where  the  constants  a,b  appearing  (19)  as  well  as  in  (18)  are  in- 
troduced to  allow  us  to  take  f(0)  =g(0)  = 1,  f and  g are  arbi- 
trary functions  to  be  determined  next. 

2 

ii)  Coefficents  of  e . 

The  equation  for  takes  the  form 

9a.  99  ■)’)'> 

4£(T?T  + ak  Tn  f,)  ” 2Ria  kf'  + 3a  kff"  - b £gg'  - abifg' 


q-  t \ 

i 


- abkf'g  + 2£(ak3f"’  + 2bS-3g"'  ) = 0.  (20) 

The  terms  appearing  in  (20)  can  be  grouped  into  three  kinds: 

a)  Secular  terms:  Those  which  are  independent  of  n.  There 

h2k2 

are  three  of  them  in  (20);  ka(-2R^af'  + 3aff'  + — ^ — f'"  ). 

Upon  integrating  these  with  respect  to  n,  we  obtain  secular  be- 
havior, i.e.  a contribution  which  becomes  unbounded  in  time  or 
space.  We  set  these  secular  terms  equal  to  zero  and  obtain  after 


letting 


R1  = I ' hok2  = 3a' 


an  equation  which  f has  to  satisfy,  i.e. 

f’"  + 3ff’  - f’  = 0,  or  f = S(£)  = sech2(|) 
Similarly  we  have  from  the  equation  for  8 , 

L1  = !•  h2£2  = 3b,  g = S ( n)  = sech2 (^) . 


I 


10 

b)  Non  Local  terms. 

These  are  not  secular  in  the  present,  they  will  be  if  left 
as  they  are.  The  solution  of  due  to  them  comes  under  an  in- 
tegral sign.  Physically  this  represents  the  memory  of  a collison. 

In  (20) , we  identify  the  following  two  terms  to  be  in  this  category. 


4ka£  f abkgf  ' . 

dri 

We  again  set  this  equal  to  zero,  and  solve  for  0Q  as 


eo  = 

44  J 

-n 

g (n  x)  dn  x • 

(24) 

Similarly  we  have 

4>  = 

o 

a | 
4k 

£ 

f(51)dC  . 

(25) 

+rt> 


(The  above  result  agrees  with  an  earlier  work  by  Oikawand  Yajima.) 

The  choice  of  the  Lower  Limits  of  integrations  will  be  explain- 
ed in  the  next  paragraph.  Note  that  the  term  abkgf'  as  mentioned 

earlier  is  non  secular  at  this  order.  Had  we  let  0 = 0 we  could 

o 

have  solved  with  a term  of  the  form  f'/n  g(o^)dn^.  Such  an  ot^ 
will  cause  a secular  solution  in  o^. 

We  now  identify  the  right-  and  left-going  waves  as  those 
described  by  f and  g respectively.  We  shall  specify  the  location 
of  these  waves  before  and  after  collison  in  terms  of  asymptotic 
values  of  £ and  n as  follows: 

Before  After 

collision  collison 

Right-going  wave  f n -*■  n +°° 

Left-going  wave  g £-*■+<*>  £ • 

We  can  then  calculate  from  (24)  and  (25)  the  shifts  of  the  waves 
after  collision.  These  phase  shifts  make  a uniform  translation 


I 


I 

f 


f 


> 
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of  the  5 and  n while  leaving  the  profiles  intact.  We  shall  call 
these  simple  phase  shifts. 

c)  Local  terms. 

The  remaining  terms  in  (20)  can  be  integrated  to  give 

2 

a1(5fn)  = ■§■  b2g2  + fg  ~ Y~  g + a2F1(C)  (26) 

similarly 

2 

61(C,n)  = | a2f2  + fg  _ a_  f + b2Gi(n)  (27) 

where  F1(C)  and  G 1 (n ) are  two  arbitrary  functions  to  be  determined 
from  the  consideration  of  the  secular  terms  in  the  next  order  of 
approximation. 


iii)  Coefficients  of  e ^ . 

The  number  of  terms  in  this  case  is  already  rather  formidable, 
we  shall  list  them  separately  in  the  different  groups  only. 

a)  Secular  terms:  Setting  these  equal  to  zero,  we  have 

F£  + (3f-l)F1  = (2R2  - i|)f  + | f2  + j f3.  (28) 

The  first  term  on  the  right  of  (30)  gives  a solution  of  Fx 
as  (2R2  - i|)  (f  + -|  f')  which  is  unbounded  as  (+  +°°.  We  thus 
require  the  coefficient  of  this  term  to  vanish  i.e. 


(29) 


This  gives  a second  order  correction  for  the  wave  speed. 

The  remaining  equation  can  be  readily  solved  and  we  have 

F1  = - | f2  + f.  (30) 


: 

I 

- 

i 


. 

i 
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We  have  now  the  complete  solution  for  and  3^: 

,,  \ 1/71,2  2 2,2.  1 ,,  2 , 2,.  ab  , . 

a1(C»n)  = g-(7b  g - a f ) - ^(b  g - 2a  f)  + j-  fg,  (31) 

n if  \ I,,  2,2  ,2  2.  1,  2,  ,,2  . ab  , /00. 

62(5,11)  = -g  (7a  f.  - b g ) - a f - 2b  g)  + y fg.  (32) 

To  the  result  of  Eq.  (30)  we  can  also  add  one  of  the  bounded  homo- 
geneous solutions  of  (28),  which  is  the  first  derivative  of  f. 
However  as  we  go  to  higher  order,  it  is  seen  that  such  an  addi- 
tion only  amounts  to  a uniform  shift  of  the  origin  of  5 which 

represents  a simple  phase  shift  defined  above.  We  will  henceforth 
drop  it  entirely, 
b)  Non  Local  terms. 

These  terms  will  give  the  solution  of  0^  and  <j>  ^ as 
»i <e ,n > = |f^  g(n1)dn1  + g2-!^  + |)bgldn1(  (33) 

1«'">  - Si*  f «5i>d  1 + Si/  a2f2  - 'F  b + f>a£>«r 

(34) 

Each  of  the  second  terms  above  is  of  similar  nature  to  the  first 
order  phase  shift  formulas  (24)  and  (25).  They  give  simple  phase 
shifts.  Note  however  that  the  first  term  in  (33)  depends 
on  C as  n +°°.  Since  0^  enters  into  the  argument  £ of  the  func- 
tion f,  we  see  that  the  wave  profile  of  the  right-going  solitary 
wave  differs  from  the  one  before  collision.  It  tilts  backward 
to  the  direction  of  propagation  of  the  wave.  (See  Fig.  1 and  the 
discussion  in  ii  below) . Similar  behaviour  appears  in  the  left- 
going wave  because  of  the  first  term  in  (34) . 
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In  section  V we  will  study  the  time  evolution  of  these  un- 
symmetrical  waves  and  show  that  its  unsymmetrical  part,  acquired 
during  collision  will  separate  and  propagate  away  from  the 
symmetric  waves  of  permanent  type  which  are  identical  to  the  waves 
before  collision. 

c)  Local  terms. 

After  integrating  the  local  terms,  we  obtain 

a 2 = J2  b3g3  + 2afa2f92  ' r*2bf29  " f ab2f9  + I a2bf9  + 


||  b3g2  - b3g  + a3F2(U-  (35) 

A similar  expression  is  obtained  for  f$2,  i • e • 

62  = 32  a3f 3 + 2a2bf2g  - | ab2fg2  - | a2bfg  + f ab2fg  + 

43  3_2  7 3^  .3  r . > 

y2  a f - -jo  a f + b G2(n)  ' (36) 

Thp  arbitrary  functions  to  integration  F2(5)and  G2(n)  like  its  counter- 
part in  the  previous  order  of  approximation  must  be  determined  from 
the  secularity-f ree  condition  in  the  next  order  of  approximation. 


4 

iv)  Coefficients  of  e . 

We  shall  focus  only  on  the  secular  terms  in  this  order. 
Setting  them  equal  to  zero,  we  obtain  an  equation  for  F2  similar 
to  (28)  as  follows 


Y"2  + ( 3f-l)  F 2 


(2R3"  H)f 


393  f2 
~ 160  f 


591  ,4 
’ 34“  f 


(37) 


The  above  equation  is  solved  by 


R - 55 

3 “ 112 

and  f2<c>  - fi  £3 

- ill  f2  + ” f 

160  1 40 

(38) 

Similarly 

r - 55 

L3  ~ 112 

and  V'”  - iU  93 

217  2 ^ 43 

160  g 40  g 

(39) 

for 

the  left  going  wave.  This  completes  our  result  for 

“2 

and  6 

as 

follows 

a ~ 
2 

1 .3  3 

- jj  b g + 

2ab2f g2  - i a2bf2g  - 

| ab2f g + | a2bf g + 

43 

32 

j_3  2 
b g 

7 

- To  b 9 

43  3 217  3-2 

+ 40  a f " Ho  3 f 

+ liZ  a3f3 

160  f ' 

(40) 

and 

B2 

= i-  a3f 3 + 
32  a f + 

2a2bfg2  - j ab2fg2  - 

§ a2bfg  + | ab2f g + 

43 

32 

2 2 

- 3,  a3f 

+ il  h3„  - 217  .32  . 

197  . 3_3 

/ A 1 v 

III-B . Summary  of  Results  in  this  Section. 


We  have  obtained  in  this  section,  the  following  results 
for  the  velocity  at  the  bottom,  W,  the  perturbed  free  surface 
elevation  £,  the  wave  speeds  Cn,CT  , and  the  phase  functions  0 ,<J> 

W(f,g)  = eC{  af  - bg  + e [-a2f 2 + b2g2  + ^(a2f-b2g)  ] + 


2 , 6 , 3 -3  ,3  3.  27.  3,2  ,32,  71.  3,  ,3. 

e [^-(a  f -b  g ) - JqU  f -b  g ) + jjj(a  f-b  g) 


+ |(a2bfg-ab2fg)  - | (a2bf 2g-ab2f g2 ) ] } 


C(f,g)  - e{af+bg+  [|-(a2f 2+b2g2)  +^abfg  + -^(a2f+b2g)J  + 


2.101,  3-3,3  3,  1 , 3-2  ,3  2,  ^ 7,  ,2-  2 , 

+ e l~80(a  f +b  g ) - jg(a  f _b  g ) + ^(ab  fg  + 

+ a2bf2g)  - |(a2bfg  + ab2fg)  + |(a3f+b2g)]}  (43) 


1 + Li  + li  22  + 15_  33 
1 2 40  £ 112  £ a 


, . eb  . 19  2,  2 , 55  3,  3 

1 + 2-  + 40  £ b + 112  £ b 


-*r«- 


bg  - ^ ae  - |e)g  dn 


+ 9eaf(uJng  dn}  (46) 


<P  = §*{]  (1  + |eaf  - be  - |S-)  f d£  + 9ebg(n)[  f dU  (47) 

■f  00  4*°0 

2 2 

where  f = sech  £/2,  and  g = sech  n/2  with  5 and  n as  defined 


by  Eqs.  (10)  and  (11) . 


- — 4- 
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To  compare  these  results  with  the  earlier  works  on  one  soli- 
tary wave  of  Laitone,  Grimshaw  and  Fenton,  we  note  that  £ 
as  used  here  is  different  from  their  non  dimensional  wave  ampli- 
tude parameter,  which  is  defined  as  the  ratio  of  maximum  wave 

a 

height,  a , to  the  undisturbed  depth  (i.e.  — ^ — ).  They 

o 

expressed  the  wave  number  and  wave  speed  in  powers  of  this 

parameter.  On  the  other  hand  we  , in  this  paper,  have  fixed  the 

1/2  . . 
wave  number  parameter  as  £ and  expanded  the  other  quantities 

in  terms  of  e.  From  (44)  we  see  that 


2.2,5, 


£ (1,0)  = ea  + a e (j)  + aJe 


3.3  13 


+ 0(e) 


anu  o o c i i ii  4 

£(0,1)  = eb  + b ( j)  + bJeJ  + 0(e*)  . 

Defining  £_  = £(1,0)  and  £ = £(0,1),  each  of  which  is  the 

K Li 

amplitude  parameters  used  by  Laitone,  Grimshaw  and  Fenton,  we  have 
after  inverting  the  above  relationship 

5 2 3 3 , 4. 

“ ' ' 4‘R  + 2 ER  + °‘V  (49) 

€b  " el  • I eL  + 5 el  + 0|el>  • 

We  now  rewrite  (44-46)  in  terms  of  eR  and  eL  as  follows: 

ir  \ - , 3 , .2  2,101  .3  151  f2  5 , . 

£<S'n)  eR^f  + 4 eR^f  + eR  ^ 80  f 80  f 8 ' 

, 3 , 2 , A 2,101  3 151  2,5  . , 

eL{g  + T EL(g  -g)  + eL( gQ-  g - qq-  9 + i g,}  + 

eReL{l+  l?(eRf+eLg)  - Ii(£R+eL)1}fg  (50) 

CR  . . eR  3 2 3 3 ' 

c“  = 1 + T - 20  eR  + 56  eR  (51) 


Zk  = ! + Ik  _ 3_  2 + 3_  3 
C 1 2 20  eL  56  eL 
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~r1(1  - ! eR  + m 4){*-V  + T2  hf'1  + eL(l  5 - 1) 

O f -*-oo 


- j±eR]g  d n + 9eRf  I g dn]} 


rj  4 (1  8 eR  + 128  eR)  {x  " CRt  + 0} 


n l 
2 = IT 


>N  - I el  + m el)(*  + V * I2  T^J1  * er<!£  - 5> 


Y-  eL]f  d£  + 9eLg  j^f  de  ] } 


ho\J  4 (1  " S'  eL  + 128  eL){x  + CLt  + t}  * 

The  last  list  of  (53)  and  (54)  serve  to  define  0 and 

The  following  qualities  are  presented  for  convenience  of  com- 
parison with  the  experimental  results  of  Chan-Street,  and  Maxworthy. 


(i)  phase  chages 


Before  collision 


A0  = 


n -oo  g -*■  o 

£ •+  +oo  ^ o 


n -*•  0 = h 


1 + 9ho  3^  e Rf 


f n " 0 - hoJ^ [1+ir  - + 9hof-r  *Rf 

After  collision  / 

[5  , , .hJjiEtl+l£  - 11El1  -9hopeL9 

40  ’ ho^(1  * r ' rER>  + 9ho'j?ERf  (55» 

4,"-hov  3s'1  + T-  - I1el>  - <56> 

The  first  term  in  the  right  of  (55)  or  the  first  term  in  the 
right  of  (56)  represents  a simple  phase  shift  (a  uniform  translation 
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without  change  of  wave  profile) . This  agrees,  in  the  first  order, 
with  that  of  Oikawa  and  Yajima. 

(ii)  Distortion  of  Wave  Profile. 

The  second  term  in  the  right  of  (55)  and  (56)  is  a function  of  £ 
or  n respectively,  i.e.  different  phase  shift  at  points  in  the 
wave.  This  causes  a distortion  in  the  wave  profile.  For  waves  of 
large  separation,  the  localized  interaction  terms  which  are  pro- 
ducts of  f and  g)  in  (50)  vanish.  Following  the  right-going  wave 
we  then  have  after  setting  g = 0, 

C - eR[f  + 4 eR(f  f)  + £r(80  f 80  f + 8 f)^'  (57) 

with  f = sech2[hoJ^(l  " | er  + §8  £r)  <x  - cRt  + 9)]»  (58) 

where  the  values  of  ® before  and  after  collision  are  given  in  (i) . 

In  fig.  1 we  plot  the  perturbed  free  surface  elevation  before  and 
after  collision  according  to  Eq.  (57)  and  (58) . Before  a collision, 
0=0,  the  profile  is  symmetric  but  thinner  than  the  first  order 
approximation  given  by  (58).  After  the  collision,  the  wave  becomes 
unsymmetrical  and  tilts  backward  to  the  direction  of  its  propagation. 
We  have  also  plotted  the  difference  between  the  unsymmetric  and  the 
symmetric  profile  (graph  No.  3) . The  propagation  of  these  tiled 
waves  will  be  discussed  in  the  next  section. 

In  addition  to  the  tilting  discussed  above,  a second  hump  will 
appear  in  the  right  branch  of  the  wave  profile  for  values  of 
-1/4 

eD  > (48)  s .3799.  However  for  such  an  e the  higher  order 

terms  which  are  neglected  in  © becomes  comparable  to  those  terms 
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returned.  The  significance  of  this  hump  is  thus  not  clear  and 
we  shall  say  no  more  about  it. 

(iii)  Maximum  (Run-up)  Amplitude  during  Collision. 

For  two  head-on  colliding  solitary  waves  with  their  maximum 
heights  before  collision  defined  as  e and  e , the  Run-up  at  a 
point  (£,n)  is  defined  by  the  value  of  the  perturbed  free  surface 
elevation  of  Eq.  (50) . It  is  easy  to  see  that  the  maximum  Run-up 
exists  at  the  point  (£,n)  where  f = 1,  g = 1 and  hence 

CReL  3 

Maximum  Run-up  = e + e + — = — + e e (e_+eT ) . (59a) 

K J_»  Z oKLiKIj 

For  two  identical  solitary  waves  eR  = el  and  the  maximum  Run- 
up is 

2 

0 

Maximum  Run-up  = 2er  + ^ er  . (59b) 

In  figure  2 we  plot  the  max  Run  up  of  (59-b) , for  the  first, 
second  (Byatt-Smith) , and  third  order  approximation,  respectively. 
The  existing  numerical,  (Chan  and  Street)  as  well  as  the  experi- 
mental results,  (Chan  and  Street,  Maxworthy)  are  also  presented  (in 
the  plot) . It  is  seen  that  our  third  order  result  is  in  complete 
agreement  with  the  numerical  values  of  Chan  and  Street  for 
e _ < 0.5.  The  experimental  data  has  a sizable  spread.  Even  then, 
the  third  order  result  seems  to  represent  them  best  if  one  excludes 
those  data  of  Maxworthy  which  represent  wave-wave  interaction  ex- 
periment. 
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V.  Time  Evolution  of  the  Unsymmetrical  Waves. 

It  is  seen  from  the  results  of  the  previous  section  that  the 
after-effect  of  a collision  on  a solitary  wave  is  manifested  only 
by  the  phase  functions  9 and  4>  As  the  collision  process  comes 
to  an  end  designated  by  n -*■  + °°  and  £ -*•  for  the  right-  and 
left  going  wave  respectively,  we  see  that,  say  for  the  right-going 
wave,  all  the  terms  associated  with  g(n)  vanish,  except  those 
appearing  in  the  phase  function  wheie  g is  under  an  integral  sign. 
After  dropping  all  terms  which  are  products  of  f and  g,  we  are 
left  with  a solitary  wave  propagating  with  a constant  speed.  This 
wave,  however,  is  not  a solution  to  the  equation  for  wave  propa- 
gating without  change  of  speed  and  shape.  In  this  section  we 
consider  the  slow  time  evolution  of  this  unsymmetrical  wave. 

We  first  derive  the  appropriate  governing  equation,  and  then 
solve  this  with  the  wave  profile  emerging  from  the  collision  as 
an  initial  data. 

Since  the  £ dependence  of  0 first  occurs  in  the  third  order 
of  e,  we  take  that  a slow  time  dependence  does  not  come  in  until 
that  order  of  approximation.  We  assume,  following  the  right- 
going wave  that 

C = eS(£  ) + e2C,(C  ) + e3C3(C„,x)  + ...  (60) 

w*  = eCR[S(CQ)  + e w'2U0)  + e2W'(£Q,T)  + ...]  (61) 

where  = /e  k(x-CRt),x  in  $ 3 and  w3  allows  the  slow  time 

evolution  of  the  wave  in  a moving  frame  fixed  in  the  right-going 
wave.  We  have  also  assumed  that  the  left-going  wave  is  far  away 
and  exerts  no  more  effect  on  the  right-going  wave.  We  define 
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* 

r 


the  slow  time  variable  by 


t = e k t. 


(62) 


Substituting  (62) , (61) , and  (60)  into  the  basic  equations  (6) 

and  (7) , we  obtain  equations  for  S in  the  first  order  approxima- 

2 

tion,  and  C W!>  for  the  order  in  e . The  solutions  of  these 
eqautions  are 


, 9 5 

R1  = S = sech  . 

R2  - if'  «2  -*2  + S'  =2  “ 1 s2  + I • 


l 

J 


(63) 


These  are  equivalent  to  the  solutions  obtained  by  Laitone. 

To  the  third  order  of  approximation,  we  have  in  terms  of 


2Wi 


n 


^ + (3S-1)  + 2j 


° 9wi  ce  19  9 

aT1  d^o  " (2V  ff)s  - ¥ 5 + 

3 4 

+ 12  S - 9S\ 


(64) 


We  decompose  W-j(Cq>t)  above  into  a stationary  solution  plus  a 
transient,  i.e. 


W 


3 ( " W + v«o'T> 


(65) 


These  functions  then  satisfy 
2 

d^W 

u ^3  55  12  2 3 

+ ( 3S-1) W,  = (2R,  - ~)S  - ~ S + 12  S° 

^ 3 J 56  6 


- 9S 


(66) 


and 


2 — + — + (3S-  1)  V]  = 0 
3T  ac  9^0 

The  solutions  to  (66)  are 


W3  - I s3  - ! S2  + \ S and  R3  - 112  ' 


(67) 


(68) 


- 

*■ 


•: 

1 


* 


s 

I 


I 
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which  is  equivalent  to  the  third  order  stationary  solution  of 
Grimshaw . 

We  now  consider  the  time  evolution  of  Eq.  (67)  with  an  initial 
data  being  the  difference  between  the  wave  field  of  the  right- 
going wave  immediately  after  collision  and  that  of  (£q)  of  (68),  i.e. 
Curve  No.  3 in  Fig.  1.  It  can  be  shown  that  for  any 

Vy(?o'T)  = F^0'k)  e_1  2T  +lkCo  (69) 

is  a solution  of  (67)  provided  that  we  take 

2 2 S’  (5o} 

F (?Q,k)  = ik(k  -1)  + 2ikS(£0)  + S'(50)  + 2 k sTT1~  (?0) 

3 

and  k + k + y = 0. 

This  form  of  solution  is  first  suggested  by  Jeffrey  and 
Kakutani  (1970).  Since  F(£Q,k)  is  bounded  for  all  to  have  a 

bounded  solution  Vy(50,T)  in  we  must  take  a real  k,  which  in 
turn  requires  a real  y . Therefore  the  general  solution  for  V in 
(67)  can  be  represented  as  an  integral  over  all  real  k,  i.e. 

r "i  + i*50 

V(£0,k)  = dk  A (k)  F(£o,k)e  Z ° (71) 

Loo 

where  A(k)  is  determined  by  the  initial  data  V(£q,t=0).  The 
asymptotic  behaviour  of  V(£0,x)*  for  x 00  is  readily  obtained  by 

* Since  S'  is  a solution  of  (67)  corresponding  to  y = 0,  any 
initial  data  of  this  form  will  remain  unchanged  in  time. 

The  simple  phase  shifts  are  of  this  sort.  They  are  exceptional 
to  what  is  to  be  said  below. 
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the  method  of  stationary  phase. 


lim  V(£o,t) 

X — ►OO 


2A(k  ) F (k  ,£  ) i (k  5 - 
o o o o o 


v oT 


TT 

2 


sgn  m") 


(72) 


where  4 = - (—  + 4)  and  y = -k  (1+k^) 

2 o x 2 o o o 

As  x + ” to  have  areal  kQ,  one  must  choose  large  negative 

values  of  £ , that  is  to  say,  asymptotically  the  distortion  of 
the  ricrht-going  solitary  wave  emerging  from  the  collision  propa- 
qate  to  the  left,  in  the  opposite  direction  of  propagation  of 
the  main  wave.  The  amplitude  of  this  secondary  wavelet  as  is 
evident  from  (72)  tends  to  zero  asymptotically  due  to  dis- 
persion. 

In  conclusion  we  state  that  when  two  solitary  waves  have  a 
head-on  collision,  the  identities  of  each  colliding  wave  are 
preserved,  (i.e.  their  wave  profiles  remain  unchanged).  However 
the  process  of  collision  does  leave  two  imprints  on  the  emerging 
waves : 


1)  Simple  phase  shifts  for  each  wave. 

These  phase  shifts  as  given  by  the  first  term  in  Eq.  (55) 
and  (56)  have  the  effect  as  if  the  waves  are  retarded  during 
the  collision. 

2)  Shedding  of  secondary  waves. 

Each  solitary  wave  sheds  a secondary  wave  with  a shape  of 
curve  3 in  Fig.  1.  These  waves  propagate  with  diminishing 
amplitude  in  the  opposite  direction  of  the  main  waves. 
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Appendix 


A1 


Substituting  (18)  into  (17)  we  obtain  a power  series  in 
Here  we  recast  each  order  of  alone. 


9a 


30. 


0(e)  = 0 


9n 


= o. 


= 0 


(Al) 


9a. 


0(ez)  = 0 


2 2 2 

- 2h*l*  9 p 

4£97T+^^-3 


9n 


f - 'V’o'Bo1 


da  36 

+ 4<tk  - k Bo>  + 


,2,  2 „ 2 

a h k 9 a a 

k vr[— = r—  + (•=  a - 2R.a)a  ] = 0 

9£  3 ^,.2  2 o 1 o 


(A2) 


9£ 


The  equation  for  8^  is  obtained  from  (A2)  by  replacing  a^  by  g ^ 
a © by  8 , & by  k,  0O  by  Cq,  n by  5 , a by  b,  and  by  L^. 
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The  equation  for  32  can  be  obtained  from  (A3)  in  a similar  fashion 
as  from  (A2)  . 
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where  F*  = a^,  G*  = b2G1,  G*  = b3G2  and  F*  = a3F2>  F1  and  F2 

are  arbitrary  functions  of  integration  of  (A2)  and  (A3)  for 

and  a2>  and  G2  are  the  corresponding  arbitrary  functions  of 

integration  for  6^  and  g2-  The  symbols  [...]  in  (A4)  represent 

4 

all  the  non  secular  terms  in  the  (e  ) order. 


